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We study the charge current transmitted through the correlated quantum dot characterized by a 
finite magnitude of the Coulomb interaction \U\. At low temperatures the correlations can lead to 
a formation of the spin (for U > 0) or charge (for U < 0) Kondo states which qualitatively affect the 
transport properties. We explore the influence of charge Kondo effect on the electron pair tunneling 
introducing the auxiliary two-component model which describes fluctuations between the empty and 
doubly occupied states on QD. 



I. INTRODUCTION 

Electronic transport through the correlated quantum 
dots (QDs) has recently attracted a considerable interest 
both from a point of view of fundamental research as 
well as practical applications [l|, Q • Both these aspects 
are possible because of a large degree of flexibility for 
adjusting a coupling of QD to the external leads and a 
controlled positioning of the quantized QD levels by the 
gate voltage. 

The tunability of the QD systems enabled controlled 
exploration of the Kondo effect a phenomenon first 
observed in metals with diluted magnetic impurities 
At sufficiently low temperatures T < Tk (where Tk de- 
notes the, so called, Kondo temperature) the spin of QD 
combines into a singlet state with the spins of itiner- 
ant electrons from the leads. In consequence the QD 
spectrum develops a narrow resonance at the chemical 
potential which at low temperatures enhances the zero- 
bias conductance to the perfect unitary limit value 2e 2 /h 
Q . On a microscopic level the underlying Kondo physics 
arises from an effective antiferromagnetic interaction be- 
tween the QD and mobile electrons as envisaged by Schri- 
effer and Wolf within the perturbative canonical trans- 
formation |5|. 

Recently several authors [1, 0] have pointed out that 
molecular quantum dots affected by the bosonic degrees 
of freedom such as e.g. phonons could show up (besides 
the side-band structure) signatures of the charge Kondo 
effect. This phenomenon has been previously theoret- 
ically predicted for the heavy fermion compounds [§] 
and it might occur when the bipolaronic shift lowers the 
charging energy to a negative value U <0. The essential 
physics involved in the charge Kondo effect relies on a 
neutralization of the electron pair charge at the negative 
U center (in the present context on QD) by electrons 
(or holes) from the adjacent leads. Due to a preferred 
double occupancy of QD there is activated a mechanism 
of the electron pair tunneling which manifests the quali- 
tative features of charge Kondo effect in the differential 
conductance 0, [1] and in the thermopower • 



In what follows we propose a phenomcnological two- 
component model which allows for a simple description 
of charge tunneling under the circumstances when the 
empty |0) and doubly occupied |f |) states are degener- 
ate (which is a necesarry condition for realization of the 
charge Kondo effect 0,0,11]). We discuss some prelimi- 
nary results for the QD spectral function and the differ- 
ential conductance obtained in the symmetric case using 
approximate treatment for the on-dot correlations. 



II. SPIN VERSUS CHARGE KONDO EFFECTS 



For a description of charge tunneling through the cor- 
related quantum dot we use the Anderson model 4] 
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Operators c^ a correspond to annihilation (creation) of 
electrons in the left (3 = L or right h.s. j3 = R lead. 
The energies £k/3er = £ko- ~ Hp are measured with re- 
spect to the chemical potentials which under nonequilib- 
rium conditions are shifted by an applied bias V through 
A*L — IJ>R = eV ■ The other terms containing V£J describe 

hybridization of the QD to external leads. As usually dV* 
denote the annihilation (creation) of electron at the QD 
energy level Ed and U is the on-dot Coulomb potential. 

Let us first start by focusing on a widely studied 
case of the repulsive charging energy U > 0. Condi- 
tions necesarry for a appearance of the Kondo effect 
can be explained using a perturbative treatment for the 
hybridization term H hyb = Ek, ) 3,<7[^4^<' + h.c.]. 

The unitary transformation H — e A He~ A with the an- 
tihermitean generating operator A = A — A* where 
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nates Hhyb up to the quadratic terms [5( . Since for U > 
the double occupancy of QD is energetically expensive 
one can restrict to the subspace of singly occupied states 
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when the transformed Hamiltonian reduces to the spin 
Kondo model Q 



~ Kondo 
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The QD spin operator Sd can be conveniently expressed 
through iSj" = = d^d^, an< ^ = ^{d\d^ — d^d^) 

and similarly S£p = etc. Near the Fermi sur- 

face the effective coupling J^'^ simplifies to || J^f kp = 
SdJzd+U) ^F^^kp/?' ■ I n the regime of antiferromagnetic 
coupling Jk'fk F < and at sufficiently low temperatures 
T < Tk the magnetic moment of QD is perfectly screened 
by spins of the itinerant electrons. For typical QDs the 
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value of Kondo temperature Tk 

(where T ~ 2^^^ |Vk F /3| 2 P/3(e_F) and pp(ef) is the den- 
sity of states at the Fermi level) is of the order of hun- 
dreds mK. Appearance of such Kondo resonance pinned 
at the electrodes' chemical potentials enhances the zero 
bias conductance and this behavior has been indeed ob- 
served experimentally [3[. 

Formally for the negative U model the canoni- 
cal transformation can be done in the same way. 
Since the empty and doubly occupied QD states are 
then energetically more favorable therefore one fo- 
cuses mainly on the terms describing pair tunneling 

Ek.qAo-.o-' ( J k,q 4^-^/3-^6^' + h.c.J 0. Tara- 
phder and Coleman Q have shown that in the symmetric 
case Ed+U /2 = (here assuming also V — 0) the attractive 
U < model becomes exactly isomorphic to its repulsive 
U > counterpart ([T]) under the following particle-hole 



(p-h) transformation d. 
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Out- 



side the symmetric situation the p-h transformation still 
renders the structure of ([1]) with an additional Zeeman 
field B z = 2e d + U (see Ref. for details). The effec- 
tive Hamiltonian || operating on the relevant empty and 
doubly occupied QD states can be then expressed by [8| 
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1). Using this pseudospin representation the spin Kondo 
effect can be directly translated into the charge Kondo 
effect of the model ((4]). The latter takes place in a vicinity 
of the degeneracy point ej + U/2 ~ 0. So far some of its 
signatures in the pair tunneling through the QDs have 
been partly examined in Refs [y, 0] ■ 
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FIG. 1: Spectral function of the molecular quantum dot for 
the equilibrium situation (V = 0) in the symmetric case 
Ep a i r — 0, Ed = for several temperatures. We used 
g = r_L = Tjj assuming the wide-band limit Tp = 0.001D, 
where D is the bandwidth of conduction electrons. 



III. THE TWO-COMPONENT MODEL 

To consider the situation when the empty and doubly 
occupied states of QD are degenerate we introduce the 
following auxiliary model 

H QD = E d d\da + g (& f <Mt + d]d\ b 



h r b. 
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Whenever the electron pair happens to arrive on the d- 
QD we let it be stored on the side-coupled buffer de- 
scribed by the operators b"\ They obey the hard-core 
boson relations [l2j so that at most only one electron pair 
can be allocated on this side-coupled b-QD. Similar lat- 
tice version of this model has been proposed in the solid 
state physics to account for the bipolaron superfluidity 
in the crossover between the adiabatic and antiadiabatic 
regimes [i~2j . 

Loosely speaking the correspondence between ([5]) and 
the negative U < Hamiltonian © holds via substitu- 



tions g = U, E. 
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U. This relation can be fur- 



ther supported analytically in the Lagrangian language 
after performing the Hubbard-Stratonovich transforma- 
tion which eliminates Ud^d^d^di through the additional 

bosonic fields, here denoted by 6^) , 

We will show that the two-component QD described by 
© is able to capture the charge Kondo effect [f| known 
for the negative U Anderson model [l3| near the sym- 
metric situation. The underlying mechanism is driven by 
suppressing the quantum fluctuations between the empty 
and doubly occupied states. Within the phenomenolog- 
ical model ([5]) such fluctuations are present even in the 
case of isolated molecular quantum dot i.e. for V^p = 0. 
Out of 8 possible configurations being a product of the 
fermionic states |0 >, | f>, | J,>, | tl> and the hard-core 
bosonic ones |0), |1) two of them |j> ®|0) and |0 > <8>|1) 
get mixed due to the Andreev-like interaction. A com- 
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plete set of the eigenstates can be obtained using the 
transformation Um 



\B> = sin(» |0 > <8>|1) 
\A > = coa(f) |0 > ®|1) 



cos(p) | U> ®|0), (6) 
sm(<f) | U> ®|0) (7) 

with tan(2<^) = 2 E d -E — -- ^ or suc ^- nm it Vk/3 — » the 
Green's function C<j( w ) = {(da, d^))^ acquires the three 
pole structure [3] 
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with 7 2 = 1 + 4g 2 /(2^ d - E pair ) 2 and Z = 

h _|_ e -E d /k B T + e -(B d +B P air)/fcBT _|_ e -(2B d +B pair )/fc B T-J 
/ [l + 2e~ Ed/kBT + e -(Ed+E-)/k B T + e -(E d +E+)/k B T + 
2e -(E d +E pazr )/k B T + e -(2E d +E palr )/k B Ty The spectral 

weight Z of single particle level is gradually depleted for 
a decreasing temperature and its amount is transferred 
to the bonding £L and antibonding E + levels. 

To get some insight into the many-body physics we 
employ a simple approximation based on the following 
substitution for the self-energy (see the discusion in Sec- 
tion 12.5 of the review book [3]) 
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We shall focus on the wide-band limit when the hy- 
bridization couplings Tp(w) = 27r^ k |Vk/3| 2 <5 (w — Sup) 
can be treated as constant quantities Tp. Figure [T] shows 
the spectral function pd t<T (oj) = — \mQd{u + i0 + ) ob- 
tained for the symmetric case E pa i r — 0, Ed = when 
u 2 = 0.5 = v 2 and E + = g, E_ = —g. We notice the 
three-peak structure where the middle one is sensitive to 
temperature due to transfer of the spectral weight Z. Let 
us emphasize that this behavior is typical for the Dicke 
effect known in quantum optics where the narrow/broad 
energy features correspond to the states weakly/strongly 
coupled to the electromagnetic field and they contribute 
the subradiant /superradiant emission lines [15|] . Recently 
a similar concept of the Kondo-Dicke effect has been pro- 
posed in mesoscopic physics for a set of three vertically 
coupled quantum dots |l6| . Our proposal ([5]) formally 
belongs to the same class and some details concerning its 
relation to the Dicke effect have been already in pointed 
out in the review paper (l5| . 

For the present context it is important to emphasize 
that the spectral function Pd,a{u — 0) is temperature- 
independent. In order to clarify this property let us con- 
sider the bare Green's function ([5]) which at low energies 
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FIG. 2: The differential conductance G(V) as a function of 
bias voltage V in the symmetric case for fcsT = 0.2F 'l, Tr = 
Tl and g=Ti,. Width of the central (superradiance) peak is 
proportional to temperature whereas the zero bias value G(0) 
remains constant (temperature-independent) 0. 



lu — > simplifies to G^ k,3=0 (a->) ~ Z/uj (we refer here to 
the symmetric case Ed = 0). Using the Ansatz (fTTj) we 
obtain Qd(u>) — u +iz(T Z L +T R )/2 wrncn yields the following 
fixed value of the spectral function pd !<7 (w) ~ ^ r ^ r . 
Furthermore we notice that effective broadening of a cen- 
tral peak depends on the spectral weight Z, hence comes 
its temperature dependence which is illustrated in figure 
[TJ Obviously, in a weak interaction regime \g\ <C Tp the 
three peaks overlap with one another therefore the fine 
structure related to the charge Kondo state fades away. 



IV. TRANSPORT PROPERTIES 

External bias V applied between the electrodes in- 
duces the charge and energy currents through the in- 
terface. We expect that behavior presented in figure 
[T] is going to have an impact on the low-temperature 
transport properties. To determine the charge current 
Ip(V) = -e|(X) kff ctp a Ckpa) one has to use either the 
nonequilibrium Keldysh or Kadanoff-Baym formalism 
[l[. Standard calculations [ll[ yield the steady electron 
current I(V) — Il(V) = —Ir(V) given by the Landauer- 
type formula 

2e f°° 

I(V) = —J dw[/(w - fi L )-f(u - fi R )] T(uj), (12) 

where Pl — ^r = and f(u>) — \\. + exp{ui) / HrT] . 
Correlation effects show up within such approach through 
the transmission coefficient T(u>) = ^l+t S ct Pd,a(u). 
This treatment is valid for the wide-band limit Tp ^ D 
and it can be derived, for instance from the equation of 
motion method [l7T |. 

We have explored numerically the differential con- 
ductance G{V) = -jyI(V) for a number of tempera- 
tures. Figure [5] shows a representative result obtained 
for fcfiT = 0.2r^ where a middle peak corresponds to 
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the characteristic charge Kondo feature. Its width is pro- 
portional to temperature, whereas the zero bias conduc- 



tance G(0) = 2«_ dw 



duj 



T{oj) at low temper- 



atures approaches the asymptotic value limx^o G(0) = 

TT Tl+t r T,aPd,v(°)- In the previous study 0] of the 
charge Kondo physics such feature has been predicted for 
the zero-bias conductance as a function of the gate volt- 
age. From our study of the two-component model {5} we 
find the same behavior but this issue will be discussed 
separately. 



V. SUMMARY 

We have studied the charge transport through a cor- 
related quantum dot with an effective negative value of 
the Coulomb interaction U < 0. To account for the quan- 
tum fluctuations between the empty and doubly occupied 
states we have introduced the phenomenological model 



where QD dot is coupled to an additional electron pair 
buffer via the Andreev-type interaction. Focusing on the 
symmetric case, when the Kondo effect arises in the pseu- 
dospin charge channel Q , we have examined its influence 
on the spectral function and the differential conductance. 
For both these quantities we have found the narrow peak 
whose height is constant whereas its width is propor- 
tional to temperature. This resembles the properties of 
the Dicke effect [l5| which has been recently considered 
for a configuration composed of the three vertically cou- 
pled quantum dots [HI]. It might be of some interest to 
proceed with an analysis of here proposed model outside 
the symmetric case and adopt some more sophisticated 
selfconsistent treatment for the Andreev-like interaction 
on the molecular quantum dot. 
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